Lagrangian planetary equations in Schwarzschild space— time 
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We have developed a method to study the effects of a perturbation to the motion of a test 
point-like object in a Schwarzschild space-time. Such a method is the extension of the Lagrangian 
planetary equations of classical celestial mechanics into the framework of the full theory of gen- 
eral relativity. The method provides a natural approach to account for relativistic effects in the 
unperturbed problem in an exact way. 
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The problem of a binary system bound by gravitational interaction and undergoing the influence of an external 
force has been widely studied within the framework of classical celestial mechanics [|J-Q. Apossible approach to 
the problem, taking into account first post -Newtonian effects, has been recently proposed |5j. In this paper we 
focus our attention to a binary system one body of which is much more massive than the other one and achieve 
the planetary equations, describing the time variation of orbital elements induced by an external perturbation. This 
■ is accomplished by starting from the exact solution to the equations of motion for a test particle in Schwarzschild 
space-time 0. Physically the problem to be solved is the perturbation of the motion of a point-like object of mass m 
around a non-rotating one of mass M 3> m. Since relativistic effects have already been considered in the unperturbed 
problem, the only quantity that is assumed to be small is the strength of the external force. This offers significant 
advantage over a semiclassical approach, which would have otherwise considered relativistic terms as a perturbation. 

The paper is organized as follows. In Sec. II we provide a review of the exact solution to time-like geodesic 
equations in Schwarzschild space-time. In Sec. Ill we achieve planetary equations for a generic external force. Finally 
in Sec. IV we provide two applications to show the capability of the method. In the first case the perturbation is a 
drag force due to interstellar dust while the second case concerns with an interstellar magnetic field, in which the less 
ON ■ massive body is assumed to be charged. However the procedure provides a method that can be useful to a general 
q ' situation. 
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6JT). II. UNPERTURBED TIME-LIKE GEODESICS 



Let us consider two point-like spinless bodies whose masses are m and M; if m is negligible with respect to M, 
then the motion is given by the following Lagrangian (e.g. @,||): 

— »£■ <«> 



r 2 (d9 2 + sin 2 9 dip 2 ) (2.2) 




dr 



2 



is the Schwarzschild line element, in which we have assumed M to be at rest on the origin (hereinafter c = G = 1; 
conventions and notations as in Ref. |7]]). Motion equations are derived from Lagrangian in the usual way: 

dt{dv dx~ u ' [Z - 6) 
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where x and v are the radius vector and velocity of m respectively. Assuming vector x (and v) to be expressed in 
cartesian coordinates, the solution to motion equations can be written as: 



(cosw cosfi — sinw sinfi cosi) a; 1 
(cos a; sin f2 + sin a; cosfi cosi) x 1 



(sin u cos f2 + cos w sin 12 cos i) x , 
(sin a; sin 17 — cos a; cosfi cosi) x , 



sinw sini i 1 



cosw sini x 2 , 



(2.4) 



where w, 17 and i are the usual Euler angles, defining the rotation that connects the observation reference frame with 
the intrinsic frame of the motion. In classical celestial mechanics they are usually referred to as argument of periastron 
(the angle in orbital plane from the line of nodes to the perihelion point), longitude of the ascending node (the angle 
measured from the positive x axis of the observer to the line of nodes) and inclination of the orbit (the angle between 
the orbital plane and the x-y plane of the observer), respectively M. Besides 
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r cos ( 
r sine 



(2.5) 
(2.6) 



is the solution to the problem in the particular reference frame whose a; 3 axis is normal to the plane of the orbit. In 
the parameterization of Ref. Jq] one has: 



F(9) 



t-T 



L 



dt 



a(l - e 2 ) 
1 + e cos x ' 
2 

y/l — 6fi + 2/ie 
* d\ 



< e < 1, 



.2 2' y 2. 



M 



sin 2 A 



a(l-e 2 )' 
E f x d(j) 



k 2 = 



d X (1 - 



2M\ 



2fie 



v /a(l-e 2 )M 
L 



y/a(l - e 2 )M 
E 



y/(2n - l) 2 - Ve 2 



1 



2M 



(2.7) 
(2.8) 

(2.9) 
(2.10) 
(2.11) 

(2.12) 
(2.13) 
(2.14) 



The other elements of the orbit a, e, and T are the relativistic extension of the usual Kcplcrian parameters, which 
they reduce to in the classical limit (i — 0, namely the semimajor axis of the ellipse, the eccentricity, and the time of 
periastron passage respectively. In Eqs. ( (2.8j ) and (2.11), integration constants were chosen so that = and t = T 
for x — 0. 

Furthermore, in order to allow the orbits to be confined within r\ < r < r% H, the following inequality must hold 
true 



< 



1 



2 (3 + e) 



(2.15) 



The orbits described through parameterization ( j2.7| )— ( pTl4 ) , together with inequality ( 2.15 ), are the relativistic ana- 
logues of usual Keplerian ones, to which they reduce in the limit /j, — > 0. If condition ( 2.15 ) is not met, there do not 
exist stable orbits and the body will eventually plunge into the singularity |6). 



III. RELATIVISTIC PLANETARY EQUATIONS 



The solution considered so far holds true under the assumptions that the motion of M is negligible and no external 
force is invol ved . When either of these perturbations cannot be neglected, the problem can stil l be approached using 
Lagrangian (2.1) provided that a perturbation term Q is added to the right-hand side of Eq. ( [2. 3D ||. Namely: 
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dt{dv) dx~ Q - 

If Q may be considered to be small, one is allowed to use a perturbative approach to the problem. To this aim 
we developed a procedure that is the general relativistic analogue of the usual Lagrangian planetary equations, in a 
Schwarzschild space-time jl]-Q. We start by defining the Hamiltonian in the usual way: 

n = p—-C, P=TT71xT- ( 3 - 2 ) 



Hamilton equations are thus written as: 



dt~ dp' dt~ dx + 14 {6 - 6) 

In order to find out planetary equations one requires that the solution to above equations has the same form of 
Eqs. (2.4)— ( pT4 ); obviously this implies that the orbital parameters a, e, T, u, fi, and i vary with time. Therefore, 
setting x — x(Cj,t), where Cj (1 < j < 6) are the orbital elements a, e, T, w, ft, i, we have: 

dx dTL dp dTL 

~dt = lfp~ 1 ~dt = ( ' ' 



where now 



d/L 

P=jpm ( 3 - 5 ) 



and —rf- are defined in such a way that v = =§r = 7^; this is accomplished if: 



j=i ^ 

as develop ed in Ref. Using the definition for p of Eq. ( ft.5| ), together with Eqs. (3.4) and ( |3.6| ), Hamilton 
equations ( |3.3| ) can be rewritten introducing the so-called Lagrangian brackets : 

[0j ' 0fcJ ~ flC/flC* dC h ' dCj' [ ' 

This way the equations describing the time evolution of the orbital elements (otherwise constant) induced by the 
external perturbation are: 



fc=l 



Therefore the problem is reduced to the solution of the above linear system in the variables The convenience of 
this approach owes to the independence of Lagrangian brackets upon time explicitly, as one can see through Eqs. (ft.4|). 
The classical version of this property is extensively treated in Ref. Q . This way we can calculate them at convenient 
time, such as t = T |l]||. After tedious but straightforward computation we obtain the only non-vanishing brackets: 



mM [l-8^ + 4/i 2 (3 + e 2 )] 

2a 2 - 2^) 2 - 4^ 2 e 2 [1 - fi (3 + e 2 ) 



a ' ^Tl = "oTF ... : = == 7. ~ —> (3-9) 



e, T = v jM -, 3.10 

v/(l - 2/x) 2 - 4 M 2 e 2 [1 - /x(3 + e 2 )] 5 

[M ^-mVM fl (l-e 2 ) SU " (3.11) 

-^(3 + e 2 ) 

m / M(l-e 2 ) l-2 M (3 + e 2 ) 

a , ^ =-7 3-) (3-12) 

2 V a [1-^3 + e 2 )^ 
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[a, n] 

[e, lu] = me 
\e. ti\ 



a, oj\ cos*, 



Ma 1 - n(7 + e 2 ) 



e, u>\ cos i. 



(3.13) 
(3.14) 
(3.15) 



We notice that, except for Eq. ( |3.10 ), the leading term in the above brackets is th e cl assical one, in the limit /i — > 0. 
As for [e, T], its classical limit vanishes. We are now in a position to invert system (SL8) and obtain the time evolution 
of orbital elements: 



da 
~dt 



2a - 2m) 2 - 4 M 2 e 2 ] [1 - /i(3 + e 2 )] [1 - m(7 + e 2 )] 
to (1 - e 2 )M [(1 - 6/j,) 2 - 4^ 2 e 2 ] 



de 
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| 8V^Vl-^(3 + e 2 ) ^ 
to [(1 — 6m) 2 — 4/i 2 e 2 ] ' 

y/[(l - 2 M ) 2 - 4^ 2 e 2 ] [1 - m(3 + e 2 )] [l - 2^(3 + e 2 )] 
fj,me [(1 — 6/i) 2 — 4/i 2 e 2 ] 



^3 



y/1 - /i(3 + e 2 ) [1 - 8 M + 4^ 2 (3 + e 2 )] 
to e y/Jl [(1 — 6m) 2 — 4^i 2 e 2 ] a 



^4 



dT _ 2a y/[(l - 2/i) 2 - 4 M 2 e 2 ] [1 - m(3 + e 2 )] [l - M (7 + e 2 )] 



df 



da; 
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dft 
~~dt 
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m(l-e 2 )p [(1 - 6m) 2 - 4^ 2 e 2 ] 



J 7 ! 



^[(1 - 2m) 2 - 4M 2 e 2 ] [1 - /i(3 + gj] [1 - 2 M (3 + e 2 )] 
m e/x [(1 — 6/i) 2 — 4/i 2 e 2 ] 

8v^ 



^2 



m [(1 - 6^) 2 - 4^ 2 e 2 ] 
\J\ — m(3 + e 2 ) cot z 



y/1 - M (3 + e 2 ) [1 - 8m + 4m 2 (3 + e 2 )] 
m e y^M [(1 ^ 6m) 2 — 4M 2 e 2 ] a 



T-2 



to y^M a(l — e 2 ) 

y/l-M(3 + e 2 ) 
to y / /Ia(l — e 2 ) sini 



v /l-M(3 + e 2 )cosi y/l-M(3 + e 2 ) 



■v^a(l - e 2 ) 



,01a(l - e 2 ) 



(3.16) 

(3.17) 

(3.18) 

(3.19) 
(3.20) 
(3.21) 



For any a ss igned external perturbation, it is possible to calculate its effects on orbital elements by solving sys- 
tem ( 3.16 )-( 3.21 ); as in classical celestial mechanics the first-order time variat io n of orbital elements is accomplished 
by substituting the unperturbed values into the right-hand side of Eqs. (3.16)— (3.21). 



IV. APPLICATIONS 

In this Section we provide two applications of the proposed method. In the former we consider a dissipative effect, 
which has a classical analogue. In the latter we show an effect which is merely relativistic. 



A. Drag force 



As a first application we consider the case in which the external perturbation is drag due to dust. Let us suppose 
to be within the range of velocities for which the drag force can be written as jjj : 



-f\v\v, 



f = 



C D Sp 



(4.1) 



where Cd is the drag coefficient, S is the cross-sectional area of the body, and p is the dust density, which we assume 
to be constant. Comparing Eq. (4.1) with Eqs. (3.1) and (2.14) we achieve: 
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Q = - 



fE\v\v 
(1 



2H' 



When this force is substituted in the last of Eqs. (p^q) we get 



fE\v\ 
fE\v\ 



dxi dx 1 dx2 dS^ • — "i ■) ■> 



.F5 = ^4 COSZ, 
^6 = 0. 



0Xi _ 



dx 2 _ 



(4.2) 



(4.3) 

(4.4) 

(4.5) 
(4.6) 



The first result we obtain is that from Eqs. ( |4.5| ) and (4.6) i and keep constant even in the perturbed case. The 
main feature of the perturbation is a secular decreasing of a which results in a spiraling of m around M. Since we are 
only interested in secular variation we perform an average on Eq. (3.16) over 2n in \- F° r the classical approach to 
this problem see jjj. To provide an analytically simple example, we assume the parameter e to be small. Up to the 
first order in the expansion in power of e we get: 
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d\ 



2ag/(l-7M 



1-6 



M 



9Me 



fl-6M 



« 



!M 2 / 



(1-2*0 (i-6^y 



1 - 



6Me 



«o (l"6g 



(4.7) 



where we have replaced the orbital elements in the right hand side of Eq. ( 3.16| ) with the initial ones, up to the 
first order in the strength of Q. It is interesting to notice that when the body reaches a = 6M in the motion, the 
parameterization used for the solution breaks down [see Eq. ( |2^1 5| )] , regardless the value of parameter e, and the body 
would plunge into the singularity even if there were no dust JqJ. 

This method could also be used to determine the variation of orbital elements caused by gravitational wave emission 
by the lightest body. To this aim it suffices to find the form of the back-reaction force due to the power loss. This 
force could be obtained starting from the energy emitted 



B. External magnetic field 



As a second example we study the perturbation to the motion of an electric charge q with mass m induced by an 
external magnetic field. We consider a particular magnetic field which is constant and homogeneous at infinity and 
directed along x 3 axis. The electromagnetic tensor-solution to Maxwell equations in Schwarzschild space-time with 
the above boundary condition-is herewith reported in terms of its non- vanishing contravariant cartesian components: 



p 13 

F 23 



F 12 = 



F 31 

F 



2MB 



x 2 x 3 



2MB 



The motion of a charged particle in an external electromagnetic field is described by (?J 

dx v 



d 2 x^ u dx a dx l3 \ 

dT 2 f a/3 Aj dj . J 



dr 



(4.8) 
(4.9) 
(4.10) 

(4.11) 



The only non-vanishing mixed components of the electromagnetic tensor are: 
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2MB 



2MB 
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F\ = -B 
2MB 



2MB 



x x x 3 . 



(x 2 ) 2 



Comparing Eq. fl4.1l| ) with Eqs. < p.l\) and ( |2.14| ) we get: 

1 



F k - v j 

^ _ 2M 3 



(4.12) 
(4.13) 
(4.14) 

(4.15) 



We also assume that the magnetic field is weak en ough t o ind uce only a small perturbation to the motion. This 
assumption allows one to solve planetary equations ( 3.16 )-( |3.21 ) in a perturbative way. 

We focus our attention to the orbits lying in a plane perpendicular to x 3 axis, so we set i — 0. Moreover we set 
O = 0. Later on we will prove that orbits satisfying these assumptions do exist. With this choice F\ = F\ — 0. The 
non-vanishing components of Q read: 

11 (fW + fW), 



2M 



Q 2 = 



2M 



{F\ «i - F\ v 2 ) 



(4.16) 
(4.17) 

(4.18) 

H = 0; (4.19) 

these equations imply that Eq. ( [3.21 ) is identically satisfied and fl is a constant we can take equal to zero. Besides, 
we get 



1 

Therefore by means of the last of Eqs. and assuming i — we obtain 

2M\ dt V)> 



Fa = J- 5 = — : 



2 1 



qL MB d , .> 



E 



dt 



(r 2 ). 



(4.20) 



This way we see from Eq. ( 3.16 ) that, as opposite to the classical c ase w here JF 3 = and the coefficient of T± is also 
zero, a is not a constant of the motion. If t is replaced by x Eq. ( 3.16| ) becomes the following equation in the first 
perturbation order (keeping orbital elements constant in the right hand side): 



da 

dx 



AqB J ^Md ^ V / [l-M3 + e2)][l-M7 + e2)] d 1 
V 1 ' [(1 - 6 M ) 2 - ^ 2 e 2 ] d X r 




(4.21) 



M 3 



y/1 - m(3 + e 2 ) 



a 3 (1 - e 2 ) 3 [(1 - 6/i) 2 - 4/u 2 e 2 ] d X 



r 2 o r (1 

— + 2 M r + 4 M log 



2M1 



(l-^f) 



'•() 



From the above expression it is straightforward there is no secular variation of a : in fact an average of the above 
equation over 2n in x is zero since r is a 27r-periodic function of x [see Eq. (2/7)]. Last equation can be solved to 
obtain the explicit dependence of a upon X' 



a (x) = a o + 



AeqB a 3 M y/[l - fi{3 + e 2 )] [1 - /x(7 + e 2 



[(1 - 6a/) 2 - 4/i 2 e 2 



(cos x ~ cos xo) 



(4.22) 



8qB 



M 3 



yj\ -^(3 + e 2 ) feag(l-e 2 ) 2 (cosx-cosxo) 



ag (1 - e 2 ) 3 [(1 - 6^) 2 - 4/i 2 e 2 ] (. (1 + e cosx) 2 (1 + e cosxo) 



1 + - (cosx + cosxo) 



2 e a M (1 — e 2 ) (cos x — cos xo) 



+ 4M 2 log 



(1 + e cos x) (1 + e cos xo) 
(1 + e cosx) K (1 -e 2 )-2M(l + e cosxo)] 1 
(1 + e cosxo) [«o (1 - e 2 ) - 2 Af (1 + e cosx)] 



where ao = a(xo)- The above expression together with Eq. ( [2.11 ) allows one to evaluate a(t). As one can check 
through Eqs. ( [3.17 ), (4.2C), and (4.18), e does not have secular terms either. This implies that Eq. (4.22) holds true 
at any time and hence r ranges within an interval. 
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V. CONCLUSION AND DISCUSSION 



We have derived Lagrangian planetary equations to describe the effect of a perturbation to time-like geodesic in 
a Schwarzschild space-time. Our method provides a natural way to study the evolution of binary systems, when 
relativistic effects can not be neglected, as for instance when coalescing stage is approached. 

The results we have obtained in this paper hold true for a test particle; from the physics point of view this means 
that the effect of the mass of the orbiting particle gives rise to a perturbation term Q m that can be neglected with 
respect to the one, Q ext , arising from an external perturbation. The order of magnitude of the former is expected to 
be the classical one, that is Q m ~ m 2 /r 2 . If the condition Q m <C Q ext ceases to hold, the problem deserves a further 
investigation. 

Other than the two examples considered our method can also be useful to other kinds of perturbations such as 
oblatcness or rotation of the orbiting star, energy loss caused by gravitational wave emission, or interaction with a 
third body. 
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